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Ionization corrections in a multi-phase interstellar medium: Lessons 
from a Zabs ~ 2 sub-DLA. 
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ABSTRACT 

We present a high resolution (FWHM= 2.7 km s"'), high S/N echelle spectrum for the 
Zem = 2.26 QSO J2123-0050 and determine elemental abundances for the Zabs - 2.06 sub- 
DLA in its line of sight. This high redshift sub-DLA has a complex kinematic structure and 
harbours detections of neutral (S i, C i), singly (e.g. C ii, S ii) and multiply ionized (e.g. C iv, 
Si iv) species as well as molecular H2 and HD. The plethora of detected transitions in various 
ionization stages is indicative of a complex multi-phase structure present in this high redshift 
galaxy. We demonstrate that the ionization corrections in this sub-DLA are significant (up to 
~ 0.7 dex). For example, if no ionization correction is applied, a super-solar metallicity is 
derived ([S/H] - -1-0.36), whereas a single phase ionization correction reduces this to [S/H] - 
-0.19. The theoretical impact of a multi-phase medium is investigated through Cloudy mod- 
elling and it is found that the abundances of Si, S and Fe are always over-estimated (by up to 
0.15 dex in our experiments) if a single-phase is assumed. Therefore, although Cloudy mod- 
els improve estimates of metal column densities, the simplification of a single phase medium 
leaves a systematic error in the result, so that even ionization-corrected abundances may still 
be too high. Without ionization corrections the properties of this sub-DLA appear to require 
extreme scenarios of nucleosynthetic origins. After ionization corrections are applied the ISM 
of this galaxy appears to be similar to some of the sighthnes through the Milky Way. 
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1 INTRODUCTION 

The measurement of rest-frame ultra-violet resonance lines in 
damped Lyman alpha (DLA) systems is currently the most success- 
ful technique for determining chemical abundances at high redshift 
(e.g. Prochaska et al. 2003; Wolfe, Gawiser & Prochaska 2005). 
The largest compilations of DLA abundances include column den- 
sities for over a dozen different elements in some 200 absorbing 
galaxies (e.g. Prochaska et al. 2007b; Dessauges-Zavadsky et al. 
2009). It is usually assumed that, due to the high A'(Hi) column 
density of the absorber, ionization corrections are negligible and 
that total elemental column densities can be derived from the dom- 
inant ionization state, i.e. the lowest energy species where the ion- 
ization potential is above 13.6 eV. For the majority of elements ob- 
served in absorption, such as silicon, iron, zinc, sulphur, nickel, ti- 
tanium and chromium, this is the singly ionized species, leading to 
the approximation N(X) x N(Xn). Exceptions include oxygen and 
nitrogen whose first ionization potential is close to that of hydrogen 
and whose ionization balance is governed by charge exchange so 
that the neutral atom is the dominant species. Therefore, although 
absorption from non-dominant species is observed in DLAs (such 



as C IV, Mg I and Na i) their contributions are minor and not usually 
considered in abundance determinations. The robustness of this ap- 
proach to abundance calculations has been tested numerous times 
in the literature (Viegas 1995; Howk & Sembach 1999; Vladilo et 
al 2001; Prochaska et al. 2002a), usually through models that em- 
ploy the ionization code Cloudy (Ferland et al. 1998). Although the 
exact magnitude of the corrections depends on a variety of input pa- 
rameters, most notably the shape and normalization of the ionizing 
background (e.g. Viegas 1995; Howk & Semback 1999), in gen- 
eral the literature agrees that the majority of DLA column densities 
do not require significant ionization corrections (e.g. Vladilo et al. 
2001). However, there are a few noteworthy cases where ionization 
corrections may be important (Prochaska et al. 2002a; Prochaska et 
al. 2002b; Dessauges-Zavadsky et al. 2004; Dessauges-Zavadsky et 
al. 2006; EHison et al. 2010). 

As the column density of neutral hydrogen decreases, the self- 
shielding approximation becomes less robust as more ionizing pho- 
tons penetrate the cloud. For example, it has been shown that the 
ratio of Al iii/Al 11 (or some proxy for Al 11 if it is saturated) tends 
to increase as the A'(Hi) decreases (Vladilo et al. 2001). As studies 
of quasar absorption line systems began to push down the ^(H i) 
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scale to investigate the nature of sub-DLA^J witii 1 9.0 < log ^(H i) 
< 20.3 cm"^ (e.g. Peroux et al. 2003) it was natural to re-assess the 
need for ionization corrections (Dessauges-Zavadsky et al. 2003; 
Meiring et al. 2007). Once again, it was found that, in general, ion- 
ization corrections were small (below 0.2 dex) and observed col- 
umn densities are therefore usually converted directly into metal- 
licities (e.g. Peroux et al. 2007; Dessauges-Zavadsky et al. 2009). 
Nonetheless, for some sub-DLAs, ionization corrections appear to 
be non-negligible (e.g. Richter et al. 2005; Quast et al. 2008). 



1.1 This work: methodology and context 

In this paper we present the case of a low column density sub- 
DLA towards the QSO J2123-0050 at z^bs = 2.06 whose unusual 
properties led us to re-assess the issue of ionization corrections for 
sub-DLAs. Since the methodology of this work differs from pre- 
vious studies of ionization corrections in sub-DLAs, it is useful 
to summarise our approach at the outset. The superb data quality, 
in terms of both resolution {R ~ 100,000) and S/N (up to 50 per 
pixel), has permitted a rare chance to demonstrate the presence of a 
multi-phase medium, with contributions from both a mostly neutral 
component, and a more highly ionized component. Although many 
works on sub-DLAs have previously investigated ionization correc- 
tions in sub-DLAs (see above references), these investigations have 
uniformly assumed that the gas measured in absorption originates 
from a single phase. Since the sightline towards J2123-0050 firmly 
establishes that this is not always the case, our main objective is to 
answer the following question: 'If multi-phase media are common 
in (sub)-DLAs, what will be the effect of calculating Cloudy ion- 
ization corrections with only a single phase?'. 

In order to address this question we create a range of mod- 
els with two-phase media, in which we vary the fractions of the 
'cold' and 'warm' gas (as characterised by different ionization pa- 
rameters). These models yield column densities of species such as 
Sill, Fell, AlII and AlIII. The key to our model philosophy is that 
we take these column densities as 'observed' values and adopt the 
usual empirical strategies of observers, and the assumption that the 
medium is a single phase, to attempt to derive the input parameters. 
In essence, we are attempting to recover the input model, but using 
an incorrect assumption about the ionization structure. Ultimately, 
we will quantify how wrong the derived elemental abundances will 
be under the assumption of a single phase ISM. The utility of this 
approach is that the relative contributions of components of real 
multi-phase absorbers can rarely be constrained. Multi-phase mod- 
els are therefore not usually possible in practice. Our models there- 
fore provide an indication of the likely error associated with the 
(necessary) single-phase approach. Additional uncertainties, such 
as those associated with geometry or in the atomic data, are not 
considered in this investigation. 

The paper is laid out as follows. In Section|2]we describe the 
observations and data reduction of the QSO J2123-0050 and de- 
scribe the measurement of the A'(Hi) and metal species column 
densities. In Section[3]we discuss the puzzling nature of this sub- 
DLA in the context of its ISM properties, if no ionization correc- 
tions are applied. The general fidelity of single phase models in the 
multi-phase case is quantified in Section[5]and evidence for a multi- 
phase medium in J2123-0050 discussed in Section |4] Returning 



' Different naming conventions have emerged for absorbers that exhibit 
damping wings, yet do not qualify as DLAs. The most common alternative 
name is super-Lyman limit system, e.g. O'Meara et al. (2007). 



to the specific case of J2123-0050, photoionization corrections are 
calculated and applied in Section|6] Results are summarised in Sec- 
tion|7] The molecular content (H2 and HD) of this absorber are in- 
vestigated in 2 companion papers (Malec et al. 2010; Tumlinson et 
al. in preparation). 



2 OBSERVATIONS AND DATA REDUCTION 

The sub-DLA towards J2123-0050 was observed with the High 
Resolution Echelle Spectrograph (HIRES) as part of a program 
to follow-up metal-strong absorbers identified in the Sloan Digi- 
tal Sky Survey (SDSS) by Herbert-Fort et al. (2006). These ab- 
sorbers are selected in an automated fashion which searches for 
absorption in 14 prominent resonance lines. From the ~ 20000 
quasars searched, more than 2000 metal absorption systems have 
been identified. In > 95% of cases where Lyo- is covered in the 
SDSS spectrum, the Hi column density is found to be large and the 
absorber would be classed as a DLA (Kaplan et al. 2010). Follow- 
up of a sub-sample of these metal-line selected DLAs with HIRES 
has shown that the metallicity of these DLAs is approaching the 
solar value, even at z ~ 2 (Kaplan et al. 2010). The main results 
of our high resolution survey will be presented elsewhere. In this 
paper, we focus on the case of the z^bs ~ 2.06 absorber towards 
J2I23-0050 (Zem = 2.26, r = 16.44). 

A preliminary HIRES exposure of 2 x 2700 seconds of 
J2I23-0050 with the CI decker (equivalent to a 0.86 arcsec wide 
slit) obtained on August 18th 2006, revealed the detection of H2 
and fine structure lines of carbon at z = 2.05930. Given the good 
observing conditions and the brightness of the QSO, further HIRES 
observations were obtained on April 19 2006 with a total exposure 
time of 10, 800 s but with the EI decker (0.4 arcsec wide), which 
yields a resolution of R ~ 100,000 (FWHM- 3 km s"'). Such 
high resolution is beneficial when studying the coldest phases of 
the ISM, allowing us to potentially resolve even very narrow com- 
ponents of the diffuse gas (e.g. Narayanan et al. 2006). All observa- 
tions were conducted with the blue cross disperser, an echelle angle 
of 0°, and cross-disperser angle of 1.0275°, yielding a total wave- 
length coverage of approximately 3000 - 6000 A. For calibration 
purposes, a set of standard trace flats were obtained, as well as the 
spectra of ThAr lamps (arcs) using the same instrument settings. 
We also obtained a set of pixel flats (lamp flats) at the beginning of 
the observing run to determine the pixel-to-pixel variation across 
the detectors. 

The data were reduced using the HIRedux routine which is a 
part of the XIDL packagqS The reduction involved the following 
procedures: 

(i) A flat field frame was produced from a stack of approxi- 
mately 30 flat field exposures for each of the detector's chips. 

(ii) In a similar manner, a combined trace flat frame is produced 
as a median over the series of standard flat images taken during the 
observing night. This frame is then used to define (trace) the echelle 
order boundaries (and find the order of curvature), and to determine 
the slit profile, which is used to correct the illumination pattern of 
the science frames. 

(iii) A wavelength solution is derived from the spectra of a ThAr 
lamp taken with the same setup as in the science frames. HIRedux 
performs a ID wavelength solution by fitting a low-order Legendre 
polynomial to the pixel values in the ThAr exposures versus the 



XIDL is publicly available at'http://wwvif.ucolick.org/~xavier/IDL/index.html| 
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laboratory wavelengths along the spatial centres of each order. The 
code then performs a 2D fit to all the lines from the ID solution. 
Finally, the pipeline derives the 2D wavelength map giving both 
the wavelength solution, and the line tilts for all orders over the full 
echelle footprint. 

(iv) After the raw science frame images are flattened and cos- 
mic rays are flagged, the final step in the reduction process is the 
extraction of the object and the sky spectra. The sky background is 
estimated from the pixels that fall well beyond the object aperture, 
taking into account diffuse scattered light, which is estimated by 
interpolating the pixel counts in the gaps between the pixel orders. 
After this, the procedure derives the spatial profile of the object 
point spread function, and performs an optimal extraction based on 
the order trace. 

(v) The individual exposures were coadded by combining each 
order separately, yielding a final, unnormalized 2D spectrum. 

(vi) The continuum is fitted manually using the XIDL routine 
xjcontinuum. The routine allows the user to select the parts of the 
continuum unaffected by absorption and then performs a minimum 
X^ fit on the selected data points using a spline function of a given 
order (for the HIRES data presented here the usual value of the 
spline order is around 8). 

(vii) The echelle orders are combined after normalization into a 
ID spectrum. In regions of order overlap, the orders are averaged, 
weighting by the square of median signal-to-noise ratio. The final 
spectrum has a S/N of ~ 15 per pixel at 3100 A, ~ 30 at 3500 A, 
and -40 at 5100 A. 




Velocity (km s } 



Figure 1. Fit to the Lya line with total logA'(Hi)=19.25±0.2 cm"^ (red 
line). The blue dotted line repre.sents the continuum fit and green lines are 
the 3(T bounds respectively. A two component fit (shown by tick marks) 
is required to adequately fit the asymmetric profile. The column densi- 
ties of separate components are logA'(Hi z = 2.05930)=19.18±0.15 cm"^, 
and logA'(Hi z = 2.05684)= 18.40±0.30 cm"-. The lower solid (cyan) line 
shows the 1 o" error array. 



2.1 A'(H i) determination 

To perform the column density measurement of the Lyce absorp- 
tion, the continuum and the line profile were simultaneously fit in 
a blaze-corrected section of the spectrum. The fit was performed 
using the xjitdla routine of XIDL. The Lyor profile is clearly 
asymmetric, with the red wing showing stronger damping. The as- 
symetry indicated that a two-component fit was necessary. The red- 
shifts of the components were fixed at the redshifts of the strongest 
metal line absorption in the two kinematically distinct metal line 
complexes detected in the system (see the next Section). The fit 
is presented in Figure [T] The column densities of the separate 
components are logMHi z = 2.05930)=19.18±0.15 cm^^, and 
logMHi z = 2.05684)=18.40±0.30 cm^^, yielding a total column 
density of neutral hydrogen of logA'(Hi)= 19.25 + 0.2 cm"-. Note 
that these uncertainties are dominated by systematic (e.g. contin- 
uum fitting), not statistical error. 



2.2 Metal column density measurements 

In order to derive the total column density of multi-component 
metal line complexes VPFIT 9.4j was used on the normalized data. 
VPFIT is a multiple Voigt profile fitting code that calculates a max- 
imum likelihood fitting function to the data. The code is adapted 
to fit multiple lines simultaneously, which allows for efficient iden- 
tification of blends. The goodness of fit is assessed in VPFIT by 
X~ statistics. The error estimates of the fitting parameters also in- 
clude the uncertainties induced by self-blends, as well as blends 
due to unidentified lines. The errors on individual components are 
fairly poorly constrained, but the total column density error can 
be more accurately quantified and may be quite small, especially 
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Figure 2. Fit to neutral carbon and sulphur lines towards J2123-0050 on 
a velocity scale relative to z = 2.05930. The ticks mark the position of the 
velocity components which give rise to molecular absorption. The data is 
presented in black, the error array is given in blue, and the fit is the red line. 
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Table 1. Neutral carbon and silicon column densities from Voigt profile fits. 



z 


b 

(kms-i) 


log N(CI) 
(cm-2) 


log N(CI*) 
(cm-2) 


log N(CI**) 
(cm-2) 


log N(SI) 
(cm-2) 


2.05930 
2.05955 


1.8±0.3 
7.3±4.0 


13.73±0.02 
12.91±0.03 


13.42±0.03 
12.74±0.10 


12.43±0.02 
12.16±0.11 


12.08±0.05 


Total 




13.79±0.02 


13.50±0.03 


12.62±0.03 


12.08±0.O5 



-200 -100 100 

Velocity (km s ') 

Figure 3. Metal lines towards J2123-0050 on a velocity scale relative to z = 
2.05930. The long blue ticks mark the position of the velocity components 
which give rise to molecular absorption. The short red tick marks indicate 
other fitted compoents. The data is presented in black, the error aiTay is 
given in blue, and the fit is the red line. 



if several species are fit simultaneously with the same structural 
model. Lower limits are reported for species with only saturated 
transitions and 3 cr upper limits are quoted for non-detections us- 
ing the following equations: 



W, 



ohs{3 cr) ■ 



3 X FWHM/(S/N), 



N= lA3e20 xW,/(AJxf), 



(1) 
(2) 
(3) 



where W,., and A,- are rest-frame equivalent width and wavelength, 
/ is oscillator strength, and (S/N) is signal to noise ratio at the 
observed line wavelength. Eqn[3]comes from the linear part of the 
curve of growth (e.g. Pagel 1997). 



b 

(kms-') 



logN(Sill) logN(Fell) log N(Nill) log N(AIII) log N(AIIII) log N(SII) 

(cm-2) 



(cm-2) 



(cm-2) 



(cm-2) 



(cm-2) 



(cm-2) 



log N(NI) log N(CII) log N(CII*) 
(cm-2) 



(cm-2) 



(cm-2) 



2.05684 


5.2 ± 0.2 


13.43 ±0.02 


12.93 ± 0.06 


11.86 ±0.04 








2.05698 


3.9 ±0.2 


13.20 ± 0.02 


12.62 ± 0.09 


12.05 ± 0.03 








2.05714 


4.5 ± 0.2 


13.28 ± 0.02 


12.72 ± 0.08 


11.74 ±0.05 








2.05743 


3.9 ±0.3 


13.03 ± 0.02 


12.74 ± 0.07 


10.85 ± 0.37 








2.05800 


6.2 ± 1.5 


12.46 ± 0.08 




10.81 ± 0.44 








2.05836 


8.5 ±1.3 


12.71 ±0.05 


12.14 ±0.32 


11.39 ±0.14 








2.05880 


7.3 ± 0.8 


12.90 ± 0.04 


12.75 ± 0.08 


11.65 ±0.07 








2.05896 


4.0 ± 0.2 


13.35 ± 0.02 


12.74 ± 0.06 


11.89 ±0.04 








2.05919 


5.2 ± 0.2 


13.57 ± 0.03 


12.74 ± 0.05 


12.13 ±0.03 


13.72 ±0.06 


13.38 ± 0.04 


12.55 ± 0.07 


2.05930 


4.8 ±0.1 


14.18 ±0.06 


13.18 ±0.07 


11.86 ±0.07 


14.16 ±0.04 


14.34 ± 0.02 


13.35 ± 0.03 


2.05943 


11.3 ±0.9 


14.11 ±0.06 


13.43 ± 0.07 


12.73 ± 0.02 


14.07 ± 0.07 




12.78 ±0.16 


2.05955 


6.4 ± 0.2 


13.82 ± 0.05 


13.24 ± 0.05 


12.06 ± 0.04 


14.20 ± 0.03 


14.00 ± 0.02 


13.19 ±0.03 


2.05979 


10.1 ±0.5 


13.45 ± 0.02 


13.40 ± 0.04 


11.79 ±0.06 






12.62 ± 0.06 


2.06017 


5.0 ±0.2 


13.68 ± 0.02 


13.49 ± 0.02 


11.92 ±0.04 


13.43 ±0.10 






2.06047 


6.5 ± 2.0 


12.40 ±0.11 













Total 



14.69±0.02 



14.12±0.02 



< 11.66 



> 13.51 



13.06 ± 0.04 



14.70±0.02 



14.53±0.02 



> 15.2 



13.71±0.01 



Table 2. Metal ions column densities from Voigt profile fits. The upper section (separated by double horizontal lines) shows fits for the satellite complex, the lower section for the main complex, as defined in the 
text. Totals in the bottom row are only for the main complex (lower section) 
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Table 3. CIV and SilV column densities from Voigt profile fits. The upper 
table section shows fits for the satellite complex, the lower section for the 
main complex, as defined in the text. The totals in the bottom row also refer 
only to the main complex. 



z 


b 


log N(CIV) 


log N(SilV) 




(kms-') 


(cm-2) 


(cm-2) 


2.05685 


9.0 ± 0.5 


13.32 ± 0.08 


12.70 ± 0.07 


2.05695 


6.6 ±3.2 


13.21 ±0.73 


12.07 ± 1.05 


2.05698 


4.1 ±0.4 


13.18 ±0.45 


13.00 ±0.09 


2.05704 


4.0 ±1.2 


13.15 ±0.32 


12.47 ± 0.30 


2.05711 


4.7 ± 0.9 


13.05 ±0.15 


12.68 ±0.11 


2.05721 


7.5 ±0.9 


13.04 ± 0.06 


12.34 ± 0.07 


2.05749 


19.5 ± 1.6 


13.12 ±0.03 


12.39 ± 0.04 


2.05800 


18.7 ± 1.0 


13.21 ± 0.02 


12.40 ± 0.03 


2.05838 


13.5 ±2.2 


12.60 ±0.11 


12.01 ±0.10 


2.05879 


6.3 ± 0.4 


12.75 ± 0.03 


12.22 ± 0.04 


2.05916 


55.4 ± 11.8 


13.48 ±0.17 


12.68 ±0.14 


2.05923 


6.6 ± 0.2 


13.83 ± 0.05 


13.30 ±0.04 


2.05930 


8.8 ±1.2 


13.96 ±0.18 


13.38 ±0.23 


2.05943 


28.6 ± 7.4 


13.50 ±0.10 


12.88 ±0.13 


2.05947 


12.8 ±3.8 


13.68 ± 0.34 


13.23 ±0.32 


2.05978 


12.1 ± 1.9 


12.95 ±0.16 


12.37 ±0.12 


2.06018 


6.3 ±1.2 


11.78 ±0.29 


11.83 ±0.07 


2.06033 


28.4 ± 3.2 


13.08 ± 0.07 


11.96 ±0.13 


Total 




14.48 ± 0.01 


13.86 ±0.01 



Table 4. Elemental abundances before ([X/H],-„„) and after ([X/H]r„„-) 
ionization corrections (IC(X/H). Values are for the main complex whose 
A'(Hi)=19.18±0.15 cm"-. Errors in [X/H] account for the en'or in A'(Hi) 
and N(X) and are added in quadrature. For presentation purposes, the errors 
are only quoted for the final corrected abundance, but the errors on the raw 
abundance are identical. 



The same kinematic structure (i.e. combination of Doppler 
{b) parameter and z) was used for singly ionized species, N i and 
Al m, but this model was not applicable to the more highly ionized 
species of C iv and Si iv. We also adopted an independent model for 
S I and C i, since these atoms likely trace gas in a cooler kinematic 
phase. The fits for the measured species are given in Tables [T] to [3] 
and shown in Figures[2]and[3] 

2.3 Molecular hydrogen 

Molecular transitions of hydrogen from both Lyman and Werner 
bands up to rotational level 7 = 5 are detected in two kinematically 
distinct components, at the velocity of the strongest metal complex. 
The H2 lines are located at the same redshifts as the neutral lines 
of carbon, which are also aligned with the strongest singly ionized 
metal ion lines. Fitting the H2 absorption features is complicated by 
saturation of the lowest J states and by the need for continuum and 
zero-level adjustment. Molecular HD is also detected in this sys- 
tem, only the third such detection in a DLA (see also Varshalovich 
et al 2001; Srianand et al. 2008). The HD and H2 lines are anal- 
ysed in two separate papers (Malec et al. 2010; Tumlinson et al. in 
preparation). The molecular lines are not considered further in this 
paper, but their presence is noted here for completeness. 



3 CHEMICAL ABUNDANCES, MOLECULAR 
FRACTION AND COOLING RATE 

The abundances, calculated on the basis of the raw column den- 
sity measurements (uncorrected for ionization), in the sub-DLA to- 
wards J2 123-0050 reveal a number of surprising results. These val- 
ues, calculated relative to the solar abundance pattern of Grevesse, 
Asplund, & Sauval (2007), are given in Table|4] We consider only 





log N(Sill) 
(cm-2) 


log N(Fell) 
(cm-2) 


log N(SII) 
(cm-2) 


log N(NI) 
(cm-2) 


Total 


14.69±0.02 


14.12±0.02 


14.70±0.02 


14.53±0.02 


log(X/H)o 

[X/H]„„, 

IC(X/H) 


-4.49 
+0.00 

+0.71 


-4.55 
-0.51 
+0.26 


-4.84 
+0.36 
+0.55 


-4.22 
-0.43 
-0.04 


[X/H],„,, 


-0.71 ±0.15 


-0.77 ±0.15 


-0.19 ±0.15 


-0.39 ±0.15 



the main absorption complex in this analysis, using an A'(Hi) = 
19.18. 



3.1 Metallicity Without Ionization Corrections 

Uncorrected for ionization, the metallicity of the sub-DLA (main 
complex) appears to be super-solar, i.e. [Si/H] = +0.00 and [S/H] 
= +0.36. Super-solar metallicities have previously been presented 
for a small number of sub-DLAs in the literature (e.g. Prochaska 
et al. 2006; Peroux et al. 2006; Peroux et al. 2008; Meiring et al. 
2008; Dessauges-Zavadsky et al. 2009), based on a similar analysis 
(i.e. assuming small or negligible ionization corrections). Such high 
metallicities are somewhat surprising at high redshift for a num- 
ber of reasons. First, the emission line abundances measured from 
actively star-forming galaxies at z > 2 are almost exclusively sub- 
solar (e.g. Erb et al. 2006; Maiolino et al. 2008). It is feasible that 
absorption and emission line abundances yield discrepant results 
(e.g. Ellison, Kewley & Mallen-Ornelas 2005) due to abundance 
gradients and/or local enrichment. However, although the magni- 
tude of H ii-to-H I region abundance disagreements is still debated 
(e.g. Lebouteiller et al. 2009 and references therein), it would be ex- 
pected that any discrepancy would tend towards higher abundances 
for emission lines. Moreover, emission line abundances measure 
only the metallicities of the actively star-forming regions, whereas 
absorption line measurements probe the entire galaxy along a given 
sightline and therefore yield average ISM measurements. Even in 
cases where a sightline probes regions near active star-formation, 
such as the case for GRB-absorbers, metallicities are typically 
around 1/10 Zq (Prochaska et al. 2007a). 



3.2 Relative Abundances Without Ionization Corrections 

Although the range of elements probed by our spectra is some- 
what limited, there are two notable abundance ratios of interest. 
The first are indicators of a element enhancement, which can be 
obtained from the ratios of [S/Fe]= +0.87 and [Si/Fe] = +0.51. 
The complexities of interpreting such ratios in the context of nu- 
cleosynthetic versus dust effects have been discussed by many au- 
thors (e.g. the review by Wolfe et al. 2005). What is striking in the 
case of J2 123-0050 is the extremely high values of [S,Si/Fe] com- 
pared to previous measurements (e.g. Prochaska & Wolfe 2002). 
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Figure 4. N/O ratio plotted against the system metallicity on the solar scale. 
Metallicities are determined from O/H or Si, S/H connected by the solar ra- 
tio of Si/0 or S/0 (see Table|4). The sub-DLA towards J2123-0050 (ringed 
green circle) is plotted against the DLA data (green circles) and Galactic H ii 
regions (open blue circles) assembled by Pettini et al. (2008). The dashed 
lines show predicted contributions from primary and secondary N produc- 
tion. The arrow points to the J2123-0050 ratio corrected for ionization. 
Representative error bars for HII regions and DLAs are given in the lower 
right of the figure. 



The main complication here is the unquantified effect of dust de- 
pletion, particularly for Fe ii which can raise these ratios far above 
their intrinsic values. 

A second ratio which reveals a surprising result is N/S. The 
utility of nitrogen as a cosmic clock, thanks to its primary and sec- 
ondary contributions, has been discussed extensively in the liter- 
ature (e.g. Pettini et al. 2002, Prochaska et al. 2002a; Centurion 
et al. 2003; Henry & Prochaska 2007; Petitjean et al. 2008). In 
brief, the dominant contribution of a elements to the ISM comes 
from prompt enrichment by massive stars which result in Type II 
SN. Secondary nitrogen is released on a longer timescale, once the 
seed nuclei of carbon have been established from previous genera- 
tions of stars. Hence, the secondary component of nitrogen steadily 
builds with time once the metallicity is approaching the solar value. 
The source of primary nitrogen is still a contentious issue, but the 
substantial scatter of N/a in DLAs has led to the suggestion that pri- 
mary nitrogen production occurs in low or intermediate mass stars 
and its release is hence is delayed relative to the a elements (Pettini 
et al. 2002; Henry & Prochaska 2007). 

The sub-DLA towards J2123-0050 is the first QSO absorber 
with solar or super-solar metallicity in which N has been measured. 
At such high metallicities. Galactic Hii regions show N/o- ratios 
that are correlated with metallicity (usually measured from O/H) 
as expected from the secondary production mechanism described 
above. In Figure|4]we show theratioof N/o- in the J2 123-0050 sub- 
DLA along with other DLA N measurements from the literature 
and Galactic H ii regions. The sub-DLA lies in a previously unpop- 
ulated part of the diagram - high metallicity, but relatively low N/o-. 
This is a puzzling result. On the one hand, this DLA has apparently 
experienced sufficient star formation to enrich its ISM to above so- 
lar values, but is only experiencing primary nitrogen enrichment. 
One explanation of this combination of abundances would be that 
the galaxy is both chemically young, but has experienced intense 
and productive star formation. 



3.3 Cooling rate 

The sub-DLA towards J2123-0050 is also an outlier in the cooling 
rate distribution of DLAs. Wolfe et al. (2003) derived a method to 
determine the properties of the star forming regions associated with 
the DLA sightlines through the detection of the Cii* /il335 line. 
The method estimates [CII] 158;um emission from the strength of 
the C II* /11335 absorption. Absorption resulting in the 1335 A line 
arises from ^P3/2, from which a 158 ^m photon is spontaneously 
emitted during the decay to the ^Pi/2 fine-structure state in the 
ground 2s-2p term of Cii. The [Cii] 158/jm transition is a prin- 
cipal coolant of interstellar neutral gas in the Galaxy (Wright et al. 
1991). By assuming that cooling is dominated by the 158;um line, 
the heating rate can be calculated by determining the [C n] cooling 
rate t. 

The emissivity of [C ii] 158yum can be calculated from the col- 
umn density of Cii* /1 1335 line following the expression of Pot- 
tasch, Wesselius, & van Duinen (1979): 

N(Cinhv,,A,„ 



k 



ergss Hz 



(4) 



N{HI) 

where A„; is the Einstein A coefficient, and hv„i is the energy of 
the 158yum transition. For the sub-DLA towards J2123-0050, we 
derive Z^ = 1.03 x lO^^^ ergs s"' Hz"'. Wolfe et al. (2008) have 
recently shown that l^ in DLAs follows a bimodal distribution with 
a transition at l^ ~ 10"-^ ergs s"' Hz"'. Absorbers above this crit- 
ical value exhibit higher metallicities, velocity widths and dust-to- 
gas ratios. Furthermore, Wolfe et al. (2008) suggest that galaxies 
with cooling rates higher than the critical value are located in more 
massive halos with active star formation occuring in 'bulge mode', 
i.e. removed from the gas halo. Although the sub-DLAs towards 
J2123-0050 qualitatively follows the above trend with its super- 
solar metallicity and complex velocity structure, it has a cooling 
rate that exceeds all the values in the Wolfe et al. (2008) DLA sam- 
ple by a factor of ten (although saturation of the C ii leads to lower 
limits for /^., especially at high /V(Hi)), and is also larger than the 
average Galactic disk value. 



4 EVIDENCE FOR A MULTI-PHASE MEDIUM 

The neutral inter-stellar medium is generally divided into two com- 
ponents - a warm neutral medium (WNM) with temperatures of 
several thousand degrees K and densities ~ 0.1 atoms cm"'', and 
the cold neutral medium (CNM) with T ~ 100 K and n ~ \Q atoms 
cm"^ (Field et al. 1969). Spin temperatures of high redshift DLAs 
are generally high (e.g. Kanekar et al. 2006 and references therein), 
with only one DLA exhibiting a value below 350 K at Zabs > 1 
(York et al. 2007). These high values are consistent with sightlines 
that intersect both CNM and WNM gas. Kanekar, Ghosh & Chen- 
galur (2001) show this temperature segregation explicitly for one 
DLA at z-iha ~ 0.2. At higher redshift, absorption from a range of 
high and low ionization species indicates that a multi-phase struc- 
ture in DLAs is common (e.g. Wolfe & Prochacka 2000; Fox et 
al. 2007; Lehner et al. 2008; Quast et al. 2008). In this section, we 
examine the evidence for a multi-phase medium in the sub-DLA 
towards J2123-0050. 

The sub-DLA towards J2123-0050 is a kinematically com- 
plex system with metal absorption spreading over 400 km s"' in 
velocity space requiring 15 components for a reasonable Voigt pro- 
file fit, see Figure [3] Although such a wide velocity spread is not 
unique amongst quasar absorbers, this sub-DLA is unusual in the 
range of species that are detected. In addition to the often-observed 
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singly ionized species, sucii as Siii and Fen, thie much less com- 
mon C I line is present and a rare detection of S i (see also Quasi 
et al. 2008). Substantial column densities of more highly ionized 
species are also measured, for example, for Al iii, C iv and Si iv. 
The singly and multiply ionized metal species are clearly separated 
into a main complex centred at z-jos = 2.05934 and a satellite com- 
plex at Zabs = 2.05684 (i.e. separated by ~ 250 km s"'). The neutral 
species, such as C i and S i are detected in only two components of 
the main complex. The same is true for the molecular absorption 
(Malec et al. 2010; Tumlinson et al. in prep) which is observed in 
only one (HD) or two (H2) velocity components associated with 
the main complex. The simultaneous presence of such a variety of 
ionized species for a given element and the kinematic diversity is 
strongly suggestive of a multi-phase medium in which neutral and 
molecular species occupy only a fraction of the interstellar volume. 
Further evidence for a multi-phase medium comes from the 
details of the Voigt profile fit. The neutral carbon and sulphur lines 
require a very small fc-parameter of only 1.8 km s"'. These very 
narrow lines are unresolved even in this high resolution spectra 
(R ~ 100, 000). However, because of a high range of gf values 
in which Ci transitions are observed, the b-parameters of these 
lines are constrained well. The singly ionized species are more than 
twice as broad with i-parameters ~ 4.8 km s"', such difi'erences 
are common signatures in the Galactic ISM (e.g. Spitzer & Jenkins 
1975). Ignoring the contribution of turbulence in line broadening, 
for a given chemical element that is detected in multiple ionization 
stages, the ratio of Doppler parameters scales as the square root of 
the ratio of temperatures, a factor of seven in the case of neutral 
and singly ionized sulphur. The sub-DLA studied here therefore 
presents a clear case of a multi-phase medium at high redshift. 



5 IONIZATION MODELS FOR A TWO-PHASE MEDIUM 

As discussed in the Introduction, ionization corrections are usually 
assumed to be negligible for DLAs, and have also been reported to 
be < 0.2 dex in sub-DLAs. Although this value is not large, correc- 
tions of this magnitude may be sufficient to mask nucleosynthetic 
signatures (e.g. Prochaska et al. 2002a; Quast et al. 2008). Previ- 
ous attempts to model ionization corrections have assumed a single 
phase model. Our observations of J2123-0050 clearly demonstrate 
that this is not always an accurate assumption. However, single 
phase models are useful because they do not require knowledge of 
how much of the observed A'(H i) is associated with each phase. As 
discussed by Ellison et al. (2007), it is extremely difficult to sepa- 
rate DLA hya profiles into separate component^, although species 
such as O I may provide clues in this respect (e.g. Fox et al. 2007). 
To assess the impact of the single-phase assumption we use Cloudy 
(Ferland et al. 1998) to construct a multi-phase model and then test 
how accurately the column densities can be recovered under the 
single-phase assumption. 

The modelling procedure is as follows. 

(i) A model with two phases is constructed, where the 'phase' 
is defined by the ionization parameter, i.e., the ratio of H ionizing 
photons to H atoms. The 'cold' phase gas is defined as having an 
ionization parameter of logU = -5.0, while the 'warm' phase is 
set to logt/ = -3.0. The total combined neutral column density 
of the cloud in both phases is set to be A'(Hi) =10""'** cm"^. Nine 

"* At lower column densities, asymmetries in the hya wing become more 
obvious and separation becomes easier 



versions of the model were constructed with different fractions of 
the A'(H i) in the warm and cold phases ranging from completely 
cold to completely warm gas (essentially these extrema are single- 
phase clouds). 

(ii) The metallicity of the cloud is fixed at [M/H]=-0.33 in all 
models, with a solar abundance pattern as given in Cloudy version 
07.02. 

(iii) The two-phase cloud is radiated with a mix of the Haardt- 
Madau (H&M) extragalactic spectrum (Haardt & Madau 1996) and 
the average Galactic ISM spectrum of Black (1987), which primar- 
ily affects the ions with ionization potential lower than IRyd, such 
as C I, and S i. 

(iv) For each of the nine models (which sample different frac- 
tions of cold and warm gas). Cloudy outputs the column densities 
of Fe n, S 11, Si 11, Al n and Al in separately for the cold and warm 
phase. The column densities of a given species in the cold and warm 
phases are summed to give the column densities that would be ob- 
served in a real spectrum. 

The procedure described above yields nine sets of column 
densities for two-phase models. The next step is to recover the 
metallicity of the theoretical input cloud by following the steps that 
an observer would execute in trying to model the cloud with a sin- 
gle phase. This requires adopting a metallicity and an indicator of 
log U and then running Cloudy over a parameter grid until the ob- 
served values are reproduced. All grids are calculated using a mix 
of the H&M and the average Galactic ISM spectrum stopping the 
calculations when the column density of neutral hydrogen reaches 
A'(Hi) =10""'** cm"-. The details of this stage in the modelling are: 

(i) For the metallicity, we assume [M/H] = [Feii/Hi] (where the 
column densities come from the Cloudy output). In order to test 
whether this approximation of metallicity introduces a significant 
error, we also repeat the experiment with the 'true' metallicity of 
[M/H]= -0.33 (in practice, this is unknown to the observer). 

(ii) The ionization parameter of the single-phase model of each 
of the theoretical clouds is inferred from the column density ra- 
tio of Al 11/ Al m. This ratio is commonly used as an indicator of 
ionization in sub-DLAs and DLAs (Dessauges-Zavadsky et al. 
2003; Prochaska et al 2002b). One of the empirical motivations 
for using Aln/Alm is its correlation with neutral hydrogen col- 
umn density (Vladilo et al 2001, Dessauges-Zavadsky et al 2002). 
Narayanan et al (2008) speculated that this anti-correlation might 
extend also to the lower column denisty QALs, such as weak 
Mg n absorbers. However, other work has questioned the utility of 
Al n/Al m for constraining the ionization parameter (e.g. Howk & 
Sembach 1999), whose atomic data is additionally poorly known 
(e.g. Vladilo et al. 2001; Dessauges-Zavadsky 2003). However, 
since our methodology is to adopt standard observational strate- 
gies to quantify errors in derived properties, we follow the com- 
mon practice of using Al n/Al m to constrain ionization parameter. 
Figure [5] presents the inferred value of log U for each of the nine 
theoretical two-phase clouds. The left-hand panels show the results 
of the models in which the metallicity is estimated from Fe ii/H i, 
whereas the right-hand panels show the models using the correct 
metallicity. The left-hand panels are almost identical to those on the 
right, indicating that, for these models, using Fe ii/H i as an estimate 
for metallicity is a good approximation. As expected, Al 11/ Al in in- 
creases monotonically with log U. Perhaps more surprising is that 
this approach fairly faithfully recovers the ionization parameter of 
the warm phase even if its contribution to the total neutral hydrogen 
density is as small as 10%. The difference between the real and re- 
covered log U for the warm phase is only on the order of 0.10 dex. 
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Figure 5. Derivation of the ionization parameter from aluminium ionic ra- 
tios calculated for theoretical two phase clouds, assuming a single phase. 
The top panels show a monotonic relation between Alii/Alin and the in- 
ferred log U. The lower panels show the inferred ionization parameter as 
a function of the fraction of N(H i) in the warm phase. The left hand col- 
umn is for the single-phase models that assume [M/H]= -0.33, while the 
right is for models with [M/H]=[Fe ii/H :]. Even if the contribution of the 
warm phase to the total column density of H i is as small as 10% the models 
recover the log U of the warm phase to within ~0.10 dex. 



(iii) After the ionization parameter is derived for each model, the 
ionization corrections are calculated. The fractional column density 

A'(X'+) 
N(X) ' 



of an element in a given ionization state is given as /(X'^) : 
and similarly for hydrogen f(HI) = - 
corrections /C(X/H) are then given by: 



and similarly for hydrogen f(HI) = ^^^ = ^jjr- The ionization 



/C(X/H) : 



logf^-logf^l, 



which is simply: 
/C(X/H) = log 



./•(X'+) 
/(HI) 



(5) 



(6) 



These values are subtracted from the ionic abundances to obtain the 
final abundance of an element. Comparing these corrected abun- 
dances to the input values of the original two-phase model allows 
us to assess how accurately ionization corrections for a single phase 
model recover the input abundances of a truly two-phase medium. 

The results of this experiment are presented in Figure[6]for the 
case where the intrinsic metallicity is known, and Figure [7] for the 
case which uses the 'observed' metallicity. The top panel presents 
the 'observed' ionic abundance for Si ii, S ii and Fe ii, the middle 
panel shows the ionization corrections, and the bottom panel gives 
the corrected metallicities. The black line indicates the input metal- 
licity of [M/H]= -0.33. Only the results for Sin, S n and Fen are 
shown, as these are the ions meausred in our study. For complete- 
ness, we also show in the Appendix a more complete set of models 
for a further six elements commonly measured in DLAs. The re- 
covered metallicity for all of the two-phase models over-estimates 
the input value. However, the deviation from the [M/H]= -0.33 
value is less than 0. 15 dex for both metallicity models, with the er- 
ror for Fe consistently lower than that of Si, and S. The accuracy 
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Figure 6. Ionization corrections for the theoretically modeled medium with 
[M/H]=-0.33. The ionic abundances 'measured' from the Cloudy mod- 
els are presented in the top panel for Sin (green line), Sn (blue line) and 
Fen (red dotted line). The ionization corrections are given in the middle 
panel, and the corrected abundances are shown in the bottom panel. The 
black horizontal Hne represents the input metallicity of [M/H]= -0.33. 
The single-phase based corrections recover the metallicity of the clouds to 
within ~0. 15 dex. 



of the corrections increases with increasing warm-phase fraction. 
The input metallicity is recovered when the model is either 100% 
warm or 100% cold phase, since in these cases, the single-phase 
approximation is obviously an accurate one. From these models, 
it is possible to quantify the approximate accuracy to within which 
ionization corrections can be calculated from a single phase Cloudy 
model. The experiment also shows that if even a small fraction of 
hydrogen column density is present in the warm phase, applying no 
corrections to the elemental abundances can lead to over-estimates 
up to 0.5 dex (as is the case for Si in this example). 



6 CORRECTED COLUMN DENSITIES FOR J2123-0050 

In this section we derive photo-ionization corrections for the sub- 
DLA towards J2123-0050 based on a single-phase Cloudy model 
and present the corrected elemental abundances. 



6.1 Cloudy model for J2123-0050 

As shown in Section |5] for a multi-phase medium, the corrections 
derived from the single-phase Cloudy model approximate the val- 
ues of the warm medium and become increasingly inaccurate as the 
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Figure 7. Same as in Figure [6] but for an assumed metallicity of 
[M/H]=[Fen/Hi]. 



medium becomes dominated by a cooler phase. It is unfortunately 
not possible to model the corrections for J2123-0050 with a multi- 
phase medium since the distribution of neutral hydrogen between 
the phases is impossible to constrain. For the case simulated above, 
we have shown that although ionization corrections that assume a 
single phase medium may still tend to over-estimate the intrinsic 
abundances, this effect is typically at the <0. 1 dex level. Apply- 
ing no ionization correction can lead to abundance over-estimates 
of more than 0.5 dex. We therefore take the enforced approach of 
adopting a single-phase model, but with the knowledge that fur- 
ther corrections for the multi-phase structure are likely to be small. 
Additional sources of uncertainty present in the modelling of any 
given absorber with Cloudy include assumptions for the UV back- 
ground, relative abundances, atomic data and absorber geometry. 

To calculate the ionization corrections (eqn |6} for the sub- 
DLA J2123-0050 we constrain the Cloudy model to have A'(H i) = 
10""'** cm"^, a metallicity of +0.36 (from the raw measurement of 
S n and H i), and a solar abundance pattern. In order to determine 
the ionization parameter (which, as we have shown above, approx- 
imates to that of the warm medium if multiple phases are present), 
the code matches the observed (Si ii/Al iii) ratio of 1 .63 to the values 
obtained from the model (see top panel of Figure[8}. Si n is used as 
a proxy for Al n since the only available Al ii line is saturated and 
only a lower limit to aluminium ionic ratio (Al n/Al ni> 0.45) can 
be measured from the available spectrum. 

The upper panel of Figure[8]presents the ionic ratio (Si n/Al iii) 
versus the ionization parameter from the model grid. The observed 
ratio of ions intercepts the model curve at the ionization parameter 
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Figure 8. The top panel shows the predicted Si ii to Al iii column density 
ratio from the Cloudy models with the photoionizing spectrum containing 
both average Galactic ISM and extragalactic HM spectra. The horizontal 
dashed lines is drawn at the measured ratio. The average ionization param- 
eter infen'ed by comparing the observed values to the model has a value of 
-2.46. The middle panel shows the fraction of each element in the given 
stage of ionization. Predicted ionization corrections for the selected metal 
abundances (see equation |6) from the same models are given in the lower 
panel. 



of logf/ = -2.46. The lower panels give the model inferred ion- 
ization fraction and correction curves for the ions of interest. Ta- 
ble |4]reports the ionization corrections (IC) for each element, and 
the corrected abundance on the solar scale once the ionization has 
been accounted for. The errors on the column densities are those 
determined by VPFIT (Table|2]l, and errors on abundances include 
the error in A'(H i) propogated in quadrature. No attempt has been 
made to estimate the error on the ionization correction due to the 
uncertainties mentioned above. Nitrogen (relative to Hi) is the least 
affected by ionization, which is often implicitly assumed, due to its 
ionization potential of 1.07 Ryd. Indeed, the middle panel of Figure 
[8] shows that the fraction of Ni closely tracks the Hi as a function 
of log U. Si II is the species requiring the largest ionization cor- 
rection for the limited species studied here (see the Appendix for 
corrections for other elements). 

In order to examine the influence of the photoionizing spec- 
trum and the assumed metallicity on the magnitude of the inferred 
corrections, additional Cloudy grids of varying metallicity (1/3 so- 
lar, and solar) and photoionizing radiation field (by scaling the 
Galactic ISM spectra component to 0, 1/10, 1/3, 3, and 10 of the 
original) were produced. The models with the combination of these 
parameters that acceptably reproduce the observed elemental abun- 
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dance ratios vary the magnitude of tlie resulting correction within a 
couple of tenths of a dex from the value determined by the original 
model. 



6.2 Application of ionization corrections 

We now re-examine the metallicity, cooling rate, N/o- and afFe 
abundances of the sub-DLA towards J2123-0050 using ionization 
corrected metal abundances. Although the 'raw' abundances indi- 
cated that this sub-DLA had a super-solar abundance in both Si 
and S, the corrected abundances of all elements are below the solar 
value: [S/H] = -0.19 and [Si/H] = -0.71. This sub-DLA there- 
fore remains metal-rich for its redshift, relative to DLAs at z ~ 2, 
superficially supporting the idea (e.g. Peroux et al. 2007) that sub- 
DLAs may be preferentially more metal-rich than DLAs. However, 
a recent re-appraisal of this hypothesis by Dessauges-Zavadsky et 
al. (2009) has argued that the difference in metallicity only exists 
at low redshift (Zii^s < 1.7), and have suggested that this might 
be due to a selection effect. In the case of the sub-DLA towards 
J2123-0050, the comparison with the bulk of DLAs is not a fair 
one, since this system was selected to be metal-strong. Herbert- 
Fort et al. (2006) provided the first evidence that absorbers se- 
lected based on strong metal lines have a tendency to approach solar 
metallicity. This has been confirmed by Kaplan et al. (2010) who 
find typical metallicities of metal-strong DLAs at z^bs ~ 2 in the 
range of 1/10 to 1/3 of solar. The sub-DLA studied here is therefore 
consistent with the range of metallicities measured in identically 
selected, higher H i column density systems at the same redshift. 

The determination of the cooling rate includes the A''(H i) col- 
umn density; for most Galactic interstellar sightlines the total hy- 
drogen content is indeed dominated by neutral gas so that A'(Hi) 
is a reasonable proxy in the calculations (e.g. Pottasch et al. 1979). 
However, we have shown that this is not the case for the sub-DLA 
studied here. We therefore re-calculate the cooling rates as: 



, mcinhvuiA,, _._ 



(7) 



where the total hydrogen column density determined from our ion- 
ization model is log N(H) = 20.05 cm"^. The cooling rate calculated 
for J2123-0050 decreases from /, = 1.03 x 10"-"' ergs s"' Hz"' 
to Ic = 1.39 X 10"-* ergs s"' Hz"', consistent with the 'high-cool' 
population studied by Wolfe et al. (2008). 

Finally, the corrected N/O ratio has also moved dramatically 
relative to its 'raw' values, due to both corrections in metallicity 
and N/O (determined from N/S with a correction for the S/0 solar 
abundance). The corrected value (see Figure |4j is now indicative 
of a primary plus secondary nitrogen contribution and is consistent 
with Galactic H n regions at a similar metallicity. Concern regard- 
ing ionization corrections of N/ff ratios has been previously raised, 
e.g. Izotov et al. (2001). Specifically, ionization effects have been 
appealed to in order to explain the large scatter of N/o- ratios be- 
low the primary plateau. Although we have shown that a signifi- 
cant correction is required for the sub-DLA studied here, Pettini 
et al. (2002) have argued that such corrections are unlikely to be 
responsible for all of the observed scatter. 



7 SUMMARY AND DISCUSSION 

The ISM is expected to exhibit a multi-phase structure and there 
is abundant evidence that DLAs are no exception. In this paper we 



present the case of a Zabs = 2.06 sub-DLA which exhibits a com- 
plex kinematic and ionization structure. The HIRES data obtained 
represent one of the highest resolution spectra obtained of a high z 
QSO which facilitates the study of the multi-phase ISM. Absorp- 
tion transitions detected in this sub-DLA range from species which 
trace cold gas (such as H2, HD, S i, Ci) up to highly ionized species 
such as C IV. Ignoring ionization corrections in this system leads to 
puzzling properties. 

We have shown that using a single phase model to derive 
ionization corrections in a multi-phase medium can only recover 
abundances to within ~ 0. 15 dex (for the model example investi- 
gated here). Although the exact accuracy of corrections will depend 
on a large number of factors (input ionizing spectrum, number of 
phases, their ionization parameters etc.), our experiment illustrates 
that caution is required when interpreting abundances in sub-DLAs. 
For coarse properties, e.g. distinguishing a solar metallicity galaxy 
from one a metal-poor one with Z ~ I/IOZq, the single-phase 
assumption probably does not impact greatly on the conclusions. 
However, relative abundances, which are used to identify subtle nu- 
cleosynthetic effects, may be quite susceptible to ionization errors. 
For example, although our two-phase model was constructed with 
a solar ratio of Si/Fe, the output abundances (see Figures [6] and 
|7]l yielded an over-abundance of Si by up to 0.1 dex. This might 
otherwise be interpreted as evidence for dust depletion in Fe, or 
Q--element enhancement in Si. 

In closing, we note that the sub-DLA studied here is of partic- 
ularly low A'(H i), and is unusual in the range of ionzation species 
it displays. It may therefore not be typical of the sub-DLA pop- 
ulation as a whole. However, it does alert us to the real need to 
assess the ionization state of sub-DLAs on a case by case basis. 
For example, Quast et al. (2008) recently revisited the z,ibf. = 1-15 
absorber towards HE0515-44. Originally estimated by de la Varga 
et al. (2000) to have log A^(Hi) k 20.45, Hubble Space Telescope 
spectra yield log A'(Hi) = 19.9 cm"^ (Reimers et al. 2003). As for 
the sub-DLA studied here, Quast et al. (2008) detect a range of ion- 
zation species including S i. Si i, Fe i. Si iii and Al iii which confirm 
that the sightline intersects a combination of neutral and ionized 
gas. In fact, Quast et al. (2008) conclude that the majority of the 
metal line species in the sub-DLA towards HE05 15-44, similarly 
to the one in 12123-0050, are formed in the ionized component. A 
second interesting example is the absorber at Zabs = 0.745 towards 
Q1331-I-17. Ellison et al. (2003) previously investigated the com- 
plex kinematic structure and ionization structure of this absorber. A 
HIRES spectrum obtained by one of us (JXP) also exhibits Si i and 
Fei, ions also reported by D'Odorico (2007), although neither of 
these cases has coverage of higher ionization species. The A'(H i) of 
the Zabs = 0.745 absorber towards QI331-I-I7 is not known, but its 
D-'mdex (Ellison 2006; Ellison et al. 2008) indicates that it is likely 
to be a sub-DLA. There are thus three cases where sub-DLAs ex- 
hibit absorption from neutral gas, presumably a cold phase, despite 
their relatively low A'(H i). In at least two of these cases, higher ion- 
ization species are also present, a strong indication that the sightline 
is intersecting multi-phase gas. Although it is difficult to draw ro- 
bust conclusions from such small numbers, one explanation could 
be that the ^(H i) column density of some sub-DLAs is low because 
of ionization from H i to H n, as well as a conversion of H i to H2 
(e.g. Schaye et al. 2001; Krumholz et al. 2009). Although Noter- 
daeme et al. (2008) do not find a strong correlation between A'(Hi) 
and the probability of H2 detection, the sample of the sub-DLA 
population is relatively sparse. A more extensive survey of neutral 
ions such as C i in DLAs and sub-DLAs would be most interesting 
in this regard. 
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APPENDIX A: 
CLOUDY 



EXTENDED ELEMENT STUDY IN 



Figures lAll and IA2I show the results of the Cloudy modeling de- 
scribed in Section[5]for a further six elements commonly observed 
in DLAs. These are presented here for completeness, but we note 
that these corrections are not generally applicable to the sub-DLA 
or DLA population as they are modelled for the particular parame- 
ters of J2 123-0050. 
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Figure Al. Ionization corrections for tlie theoretically modeled medium 
with [M/H]=-0.33. The ionic abundances 'measured' from the Cloudy 
models are presented in the top panel for different elemental species com- 
monly measured in DLAs and sub-DLAs. The ionization corrections are 
given in the middle panel, and the corrected abundances are shown in the 
bottom panel. The black horizontal hne represents the input metallicity of 
[M/H]= —0.33. The single-phase based corrections recover the metallicity 
of the clouds to within ~0. 15 dex. 
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Figure A2. Same as in Figure lAll but for an assumed metallicity of 
[M/H]=[Fen/Hi]. 



